Coarse-grained simulations of DNA overstretching 
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We use a recently developed coarse-grained model to simulate the overstretching of duplex DNA. Overstretch- 
ing at 23° C occurs at 77 pN in the model, about 10 pN higher than the experimental value at equivalent salt 
conditions. Furthermore, the model reproduces the temperature dependence of the overstretching force well. 
The mechanism of overstretching is always force-induced melting by unpeeling from the free ends. That we 
never see S-DNA (overstretched duplex DNA), even though there is clear experimental evidence for this mode 
of overstretching under certain conditions, suggests that S-DNA is not simply an unstacked but hydrogen- 
bonded duplex, but instead probably has a more exotic structure. 



PACS numbers: 87.14.gk,87.15.A-,87.15.La 

I. INTRODUCTION 

DNA in vivo is not just a passive molecular bearer 
of information, but is an active molecule that is able 
to respond structurally to cellular signals, be it through 
changes in solutions conditions, protein binding or the 
action of molecular machines. Some of these changes are 
mediated by very specific responses to the biochemical 
details of the protein binding, say, but in many instances 
this control is achieved through the mechanical response 
of DNA, e.g. enzymes that apply structural control by ad- 
justing the supercoiling of DNA. For these reasons, there 
has been much interest in the fundamental mechanical 
properties of DNAp^ particularly as there are now the 
means to study these properties in unprecedented detail 
using single-molecule techniques. 

One particular focus has been the response of DNA to 
tension. At low forces the B-DNA duplex responds elasti- 
cally and is well described by the worm-like chain model. 4 
However, at higher forces, typically in the 60-70 pN range 
at room temperature (although the precise value depends 
on solution conditional^), DNA undergoes a dramatic 
overstretching in which it extends by about 70% over 
a few pN. Perhaps somewhat surprisingly, even though 
the first detailed results on DNA overstretching were re- 
ported in 1996,^1 the nature of the overstretched state 
is only now beginning to be fully resolved and has been 
a controversial topic for much of this period. 

The two main proposals are that yielding corresponds 
to a transitions to a new overstretched f orm of double- 
stranded DNA (normally termed S-DNA) EE] or t force- 
induced meltingP^^ Features interpreted as pointing 
to an S-DNA mechanism include (i) the lack of com- 
plete strand separation after overstretching as might be 
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expected for force-induced melting, (ii) that the force- 
extension curves after overstretching often do not ini- 
tially follow that expected for single-stranded DNA (ss- 
DNA), (iii) that further transitions can sometimes be 
seen at higher force (that are taken to correspond to 
force-induced melting and strand dissociation) and (iv) 
the reversibility of the transition under some conditions. 
By contrast, the dep ende nce of the overstretching force 
F over on temperatur e! 7 * 13 ! and solution conditions (e.g. 
pEP and salt concentratiorP^ fits well with the force- 
induced melting picture; i.e. conditions which destabilize 
the duplex (e.g. higher temperature, lower salt) lead to 
a lower overstretching force. Furthermore, the hysteresis 
sometimes seen has been interpreted as due to the slow 
nature of reassociation after the melting of a long DNA 
molecule. 

New impetus has come to this debate from the re- 
cent fluorescence experiments of van Mameren et al. 
that clearly showed that at room temperature and at 5- 
150 mM salt overstretching occurred by unpeeling from 
the free endsFS] Although this experiment was inter- 
preted by some as proof of the force-induced melt- 
ing hypothesis (a nd h ence that the S-DNA hypothesis 
was unnecessary) P^E^l this conclusion was dispu ted by 
othersP^!!] Furthermore, subsequent experiment J 7 I 2 °H 23 
have shown this interpretation to be too simplistic, and 
provide a basis for a more balanced position. In particu- 
lar, these experiments have shown that the two modes of 
overstretching can occur depending on conditions, with 
S-DNA being favoured at low temperature, high salt, 
high G-C content and short time scales. One of the com- 
plications is that at room temperature the two mech- 
anisms can compete. So for lower salt concentrations 
unpeeling is dominant, as in the experiments of van Ma- 
maren et al\^ and the impressive high-resolution follow- 
up study of Gross et al. which resolved a saw-tooth struc- 
ture in the overstretching force-extension curves that cor- 
related with the ease with which different parts of the 
sequence could unpeel.^ However, at higher salt, an S- 
DNA mechanism becomes more feasible. For example, 
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in the elegant experiment of Paik et al. at 150 mM salt 
the transition occurred with some hysteresis when free 
ends were present (presumably due to unpeeling), but 
occurred reversibly at the same overstretching force for 
a DNA construct that had no free ends and was still tor- 
sionally unconstrained (presumably by an S-DNA mech- 
anism) . 

There has also been a lot of work using theory and sim- 
ulation that has aimed to provide further insights into 
the nature of the overstretching transition. The theo- 
retical approaches usually use simple polymer models to 
describe each of the different possible states with param- 
eters fitted to reproduce experimental behaviour! 17 * 28 ! ^ 
In particular, the wor k of Co cco and Markop^ and 
Whitelam and coworker d 17 * 34 * 35 ! has provided important 
insights into the competition between different processes, 
such as S-DNA formation versus unpeeling. However, 
as these models lack a detailed molecular representation 
of the underlying phenomena, the questions these ap- 
proaches can address is inevitably limited. 

At the other extreme, one has simulations of DNA 
stretching using fully-atomistic modelspmSB These have 
the potential to provide answers to detailed structural 
questions that are inaccessible to experiment. However, 
firstly, the computational cost of such simulations means 
that the pulling rates are often extremely fast, making 
it likely that much of what is observed is not at equi- 
librium. Secondly, it is not clear how well the atomistic 
potentials will perform under the more extreme condi- 
tions associated with overstretching. Thirdly, in some of 
the simulations, constraints or boundary conditions are 
used that restrict the range of possible behaviour, e.g. 
torsional constraints or the absence of free ends. For 
these reasons, it is perhaps not surprising that these sim- 
ulations only observe overstretching at forces well above 
the experimental values (all above 100 pN, some substan- 
tially more) and that the transitions are much broader. 
Furthermore, although they have the potential to identify 
the structure of S-DNA, no consensus has been reached. 

In the middle are simulations using coarse-grained 
models that have a simplifed representation of DNA and 
so allow longer time-scale and length-scale processes to 
be more easily accessible. These can be divided into two 
types. Firstly, simple models that attempt to obtain 
qualitative insights into the stretching behaviour that 
stem from ca pturing some of the basic physical ingre- 
dients of DNA.I2ZE2M] Secondly, there are coarse-grained 
DNA models that are fitted in order to provide a quan- 
titatively accurate description of DNA's behaviourl^MED 
For some of these models, the stretching behaviour has 
been simulated .^^] However, the observed overstretch- 
ing forces are also typically much higher than experiment 
(e.g. 320 pN in Ref . M\ 500-1100 pN in Ref.[59]and 400- 
550 pN in Ref. EH). 

Here, we simulate DNA overstretching for our recently 
developed coarse-grained DNA potential, which provides 
a quantitative description of many of the struct ural, t her- 
modynamic and mechanical properties of 

DNA L6i]62] o n e 




FIG. 1. (a) Two nucleotides represented by our model, one 
showing the rigid nucleotide unit and one the backbone and 
base regions, (b) A short DNA duplex as represented by our 
model. 

aim of the work is to provide an independent test of the 
robustness of the model by applying it to phenomena to 
which it has not been fitted. A second aim is to provide 
fundamental insights into the nature of the overstretch- 
ing transition, both by visualizing some of the mechanis- 
tic aspects of the transitions, and by deducing from the 
successes and failures of the model what ingredients are 
needed to capture different aspects of the transition. We 
proceed as follows. In the Methods section, we describe 
our DNA potential, the simulation techniques used and 
some of the theory of the thermodynamics of pulling. In 
the Results section, we focus in detail on the room tem- 
perature overstretching behaviour, before more briefly 
examining the temperature dependence of overstretch- 
ing, including looking for the re-entran t behaviour that 
has been predicted to occur near melting.^ 1 13 1 32 1 33 1 37 1 38 1 52 1 



II. METHODS 
A. DNA model 

We use the coarse- graine d model of DNA recently de- 
veloped in our groupF^SU i n particular, the model has 
been designed to provide a good description of the struc- 
tural, thermodynamic and mechanical properties of both 
double- and single-stranded DNA, all features that are 
very important if we are to be able to describe DNA 
overstretching. In addition, the wide range of appli- 
cations for which the model has been successfully used 
give us further confidence in its robustness. These 
applications so far include DNA nanotweezersj^ kiss- 
ing hairpins DNA walkerspB the formation of cruci- 
form DNA^ and the nematic transition of dense solu- 
tions of short duplexes)^ furthermore, investigations into 
hybridization, strand- exchange reactions, plectonemes, 
DNA origami and phase transitions in solutions of DNA 
junctions are ongoing. 

The model describes each nucleotide as a rigid ob- 
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ject (Fig. [lja)) that has interactions corresponding 
to backbone connectivity, excluded volume, stacking, 
hydrogen-bonding between complementary bases and 
cross-stacking. Each nucleotide has three collinear in- 
teraction sites (corresponding to the centres of the back- 
bone, stacking and hydrogen-bonding interactions) and 
a vector perpendicular to this axis to capture the pla- 
narity of the bases. We would like to emphasise that the 
attractive interactions in the model are not isotropic, but 
depend on the relative orientation of the nucleotides. It 
is these angular modulations that ensure that the DNA 
helix is right-handed and pairs in an anti-parallel man- 
ner. The full form of the potential has been given in Refs. 
I62landl64l A simulation code incorporating the potential 
is available to download.^ 

We should note a number of simplifying features of the 
model that are relevant to the current study. Firstly, we 
do not attempt to model the electrostatics explicitly, but 
instead have fitted the potential parameters for a particu- 
lar salt concentration, namely 500 mM. At this concentra- 
tion the Debye screening length is short-ranged and the 
effects of the electrostatics are included in the excluded- 
volume interactions. Fortunately, this salt concentration 
is also one of those commonly used in experiments on 
overstretching. However, as a consequence of this sim- 
plification we cannot of course examine how overstretch- 
ing depends on salt concentration. At room tempera- 
ture and 500 mM salt, experiments indicate that both the 
force-induce d melti ng and S-DNA modes of overstretch- 
ing compete! 7 * 22 ^ 23 ! 

A second of the model's simplifications is that it ig- 
nores sequence dependence in the interactions except for 
the Watson-Crick nature of the base pairing (i.e. the hy- 
drogen bonding term in the potential only occurs for 
complementary bases). Therefore, each base-pair has the 
same average interaction strength irrespective of its iden- 
tity and its neighbours. This "average-base" approxi- 
mation can be an advantage when one is interested in 
the generic behaviour of DNA, as for the most part we 
are here, but of course it excludes us from, for example, 
examining how the overstretching depends on G-C con- 
tent. A sequence-dependent version of the model is being 
developed Pi 

Thirdly, the double helix in our model (Fig. ljb)) is 
symmetrical with the helical grooves being of the same 
size. This simplification, however, is likely to be rela- 
tively unimportant for the current study. 
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FIG. 2. Schematic representation of the three pulling 
schemes used. The arrows on the DNA are in the 3'-5' di- 
rection. 



are applied to different strands, dissociation will be com- 
pletely irreversible when pulled at constant force. In our 
simulations, we pull on both 3' ends of each strand, but 
given the symmetric nature of the helix in our model and 
the simplified representation of the nucleotides, we do not 
expect any differences for pulling both 5' ends. In experi- 
ments differences between 5'-5' and 3'-3' pulling are only 
found at forces well above that for overstretching.^ 

Finally, in scheme III, as both ends of both strands are 
pulled, force-induced melting by unpeeling is suppressed. 
Therefore, if the DNA is to melt, it would have to be by 
bubble formation. Note the system in scheme III is not 
torsionally constrained, i.e. the ends are free to rotate. 



C. Thermodynamics of pulling 

Before we consider our pulling simulations, it is impor- 
tant to understand the effect of impo sing a force on the 
thermodynamics of our systems! 10 * 31 ^ The force provides 
an additional source of work being done on the system. 
Assuming the force F acts along the z direction, this gives 
an additional term Fdz to the change in internal energy 
and hence to the change in Hclmholz free energy dA, and 
A is therefore a function of z. However, it will usually 
be more convenient to consider a system under constant 
force, and so we introduce the free energy A' = A — Fz to 
describe such a constant-force ensemble. It follows that 
the change in free energy on applying a force is 



B. Pulling Schemes 

The three ways of pulling DNA that we will consider 
are shown in Fig. [2j The one we will focus on is scheme I 
where a force is imposed on both ends of the same strand, 
and the free ends allow the duplex to overstretch by un- 
peeling. Scheme II is similar in that a force is exerted on 
only one of the two strands at each end of the duplex, 
thus again allowing unpeeling. However, as the forces 



A'(F) -A'(F = 0) 



z{F)dF. 



(1) 



Therefore, the force-extension curves of dsDNA and ss- 
DNA determine the relative stabilization of these two 
forms by force. Furthermore, assuming overstretching 
occurs by force-induced melting, one just needs a correct 
description of the zero-force thermodynamics and the 
force-extension curves to correctly predict overstretching. 
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Force-induced unpeeling of DNA occurs at the force at 
which the free energy per base pair of the ssDNA and 
dsDNA states are equal, leading to coexistence of the 
two forms. In other words, at this overstretching force 
F oveT (T), the average free energy change for a duplex to 
unpeel by one more base pair is zero. We wish to em- 
phasise that overstretching, in the force-induced melting 
picture, is not determined by an equilibrium between ds- 
DNA and fully-dissociated ssDNA, but rather by when 
the duplex becomes unstable with respect to unpeeling. 
In this sense overstretching is distinct from melting, as 
the translational entropy gain from complete dissociation 
of the two strands plays no role. However, in the bulk 
limit of infinitely long chains, the contribution of trans- 
lational entropy is negligible and melting and unpeeling 
transitions are equivalent. In this limit, the melting tem- 
perature T m (F) can be defined as the temperature at 
which the average free energy change for the loss of one 
base pair from a duplex is zero at a given force F, and 
the curves F ovel (T) and T m (F) then give equivalent rep- 
resentations of the overstretching transition in the force- 
temperature plane. However, for finite chains T m (F) will 
depend on concentration, and this equivalence will be 
broken. 

We also note that care should be taken when apply- 
ing this approach to the different schemes in Fig. [2j For 
scheme I only one of the strands is force-bearing when 
single stranded. By contrast, in scheme III the force on 
each single strand is F/2. Finally, in scheme II, equi- 
librium is not well-defined, as for any non-zero force the 
most stable state is two dissociated single strands because 
they can then increase their separation without limit. 



D. Simulation methods 

We use a mixture of Monte Carlo and Brownian dy- 
namics to simulate DNA stretching for our model. We 
use Monte Carlo to obtain the equilibrium behaviour of 
the system at a particular force, and Brownian dynamics 
to study the stretching dynamics of DNA when the force 
is increased or decreased (linearly) as a function of time. 

To aid the equilibration in our Monte Carlo simula- 
tions, we use the virtual move Monte Carlo approach of 
Whitelam and coworkers.^ This algorithm is a type of 
"cluster-move" Monte Carlo that allows clusters of nu- 
cleotides to be moved at each step, where the clusters 
that are constructed reflect both the configuration and 
the proposed move. The collective motion of nucleotides 
that this algorithm introduces allows more efficient sam- 
pling for our model compared to Monte Carlo with single- 
particle moves. 

For the dynamics simulations we use the simple Brow- 
nian thermostat introduced in Ref. [7TJ The system is 
evolved according to Newtonian dynamics, but every 
"Newt time steps a fraction p of the particles have all three 
components of their velocities drawn from the Maxwell 
distribution. On time scales longer than priNowt, the dy- 
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FIG. 3. Force-extension curves for DNA at 23°C for pulling 
scheme I. In (a) results from a series of Monte Carlo simu- 
lations at constant force are presented. For comparison the 
curve for single-stranded DNA has been added. In (b) results 
from dynamics simulations are presented at a number of dif- 
ferent pulling rates. The pulling rates are given in pN/time 
step. These units can be converted into pNns -1 by multi- 
plying by 1.17 x 10 5 ns, but as with any coarse-grained model 
absolute values of time should be treated with caution. Using 
this conversion, our slowest and fastest rates correspond to 
0.0284 pNns -1 and 5.69 pNns" 1 , respectively. 



namics is Brownian. 

Modifying these simulation algorithms to incorporate 
the forces exerted on the terminal bases consistent with 
the different schemes in Fig. [2] is straightforward. To 
perform Monte Carlo at constant force it is necessary 
to include an additional term in the energy — J^i^Zj, 

i 

where Fi is the external force imposed on nucleotide i. 
In Monte Carlo, this gives rise to an extra term in the 
Boltzmann factor in the standard Metropolis acceptance 
criterion, namely exp(/3^ FiAzj). For the Brownian dy- 

i 

namics, implementing the constant force is just a matter 
of adding the external forces to the forces acting on the 
relevant nucleotides. 



III. RESULTS 

In all cases, we consider a 100 base-pair duplex with 
a random sequence. We first consider pulling by scheme 
I in detail. Force-extension curves at "room tempera- 
ture" (23° C) are presented in Fig. [3] for Monte Carlo 
simulations at different forces, and for dynamics sim- 
ulations at different pulling rates. This data shows a 
number of clear features. Firstly, the MC data (which 
we expect to be closer to "equilibrium") shows a clear 
and narrow overstretching transition at a force of about 
77 pN. At the slowest pulling rates, the dynamics results 
also show a clear overstretching plateau at a very similar 
force. However, as the pulling rate increases, there is an 
increasing tendency to overshoot this transition, and for 
overstretching to start at a higher force and for the tran- 
sition to be spread over a wider range of force, simply 
because the time scale associated with the overstretch- 
ing transition is no longer fast compared to the pulling 
rate. These effects of pulling rate are also relevant to 
the all-atom simulations as they provide an indication 
of the consequences of being increasingly far from equi- 
librium. Indeed, our results at faster pulling rates some- 
what resemble the force-extension curves seen in all-atom 
simulations P^HSU 

Secondly, beyond the overstretching transition, the 
force-extension curves follow that for ssDNA, thus indi- 
cating that DNA overstretching in our model is a result 
of force-induced melting. This conclusion is confirmed 
in Fig. [4] which depicts typical configurations as the sys- 
tem passes through the overstretching transition. From 
these shapshots, it can be clearly seen that the increase 
in extension is a result of DNA unpeeling from the free 
ends. In particular, there is no sign of any S-DNA-like 
state. Experiments suggest that the transition from B- to 
S-DNA is able to occur more rapidly than force-induced 
melting.^ 7 1 22 1 23 1 Therefore, if S-DNA exists for our model, 
one might expect that it is more likely to be seen in the 
faster pulling simulations where the overstretching is ob- 
served at higher force. However, in all the simulations 
depicted in Fig. [3jh) we never found any evidence of an 
S-DNA-like state. 

There are a number of more minor features that are ap- 
parent from the configurations in Fig. [4j Firstly, imme- 
diately beyond the overstretching transition, the duplex 
does not fully dissociate, because for short (e.g. 4 base 
pairs) intermolecular helices, the force-bearing strand can 
still approximately allign itself along the direction of the 
force, and so the extension gain associated with the loss of 
the base pairs is lower for these short duplexes. There- 
fore, a somewhat higher force is required to dissociate 
these last few bases. This effect can be seen in Fig. [5^ a) 
where overstretching leads to the loss of most, but not 
qute all, of the energy associated with intermolecular hy- 
drogen bonding — only at slightly higher forces does this 
term go to zero. 

Secondly, from a careful inspection of Fig. [4j one can 
see that the orientation of the bases in the overstretched 
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FIG. 4. Snapshots of typical DNA configurations as the 
system passes through the overstretching transition at 23° C 
for scheme I. The snapshots are from the dynamics run at a 
pulling rate of 2.43 x 1CT 7 pN/time step. 
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FIG. 5. Contributions to the potential energy from (a) 
hydrogen bonding and (b) stacking as a function of force. 
In (b) the contributions from the two individual strands is 
included. k B T/e = 0.1 at T = 300 K 



force-bearing strand is not random, but there seem to be 
short runs of stacked bases. Although completely stacked 
ssDNA in our model is helical, it is possible for short sec- 
tions of 3-4 stacked bases to orient the backbone approxi- 
mately along the direction of forceP^ For this reason, the 
force-bearing strand is able to retain a significant frac- 
tion of its stacking interactions after overstretching. Fig. 
[5jb) indicates that just after overstretching is complete, 
the stacking energy for the force-bearing strand is 57% 
of that for B-DNA and 61% of that for ssDNA at zero 
force (assuming the unpeeled non-force-bearing strand is 
typical) . 

Thirdly, it is clear from Fig. [4] and the appearance 
of intramolecular base pairs in Fig. [5ja) that the un- 
peeled strand exhibits secondary structure, in particular 
hairpins. Even though the sequence is random, hairpins 
can be formed with stems having a sufficient number of 
Watson-Crick base pairs (as well as some mismatches) to 
be stable. 
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Fig. [3^a) only provides a rough guide to the position of 
the overstretching transition, because it is hard to fully 
equilibrate the simulations in the vicinity of the transi- 
tion. One approach to locate i^ovcr, as noted in Sect. |II C| 
is to find the force at which the average free energy for un- 
peeling a base is zero. However, this approach is compli- 
cated by secondary structure formation. Firstly, it makes 
thermodynamic sampling of the transition more difficult, 
because secondary structure in the unpeeled state leads 
to free energy barriers between different states. For ex- 
ample, it might be that for intcrmolecular base pairs to 
form a hairpin must first open, or as the double helix 
unpeels the identity of the most stable hairpin in the un- 
peeled strand changes. Secondly, it obscures some of the 
basic physics of the overstretching transition, as hairpins 
will stabilize the unpeeled state, but the extent to which 
this occurs (and hence the overstretching force) will be 
sequence dependent (even for our "average base" model 
where all the interactions are assumed to be of the same 
strength) . 

Therefore, in our simulations to locate the overstretch- 
ing force, we turn off intramolecular base pairing^ 
and, but less importantly, non-native intcrmolecular base 
pairing.^] As a consequence, the free energy change for 
unpeeling every base pair is now zero (not just on aver- 
age) at the overstretching transition of our model and so 
we only need to sample a local section of the free energy 
landscape, rather than that for the whole transition. 

In Fig.[6ja) we show the free energy landscape for a 10- 
base section of the duplex for a number of forces close to 
i^ovcr- In real DNA the free energy landscape for unpeel- 
ing will have a saw-tooth structure due to the sequence- 
dependence of the thermodynamics of base-pairing (e.g. 
G-C base pairs are stronger than A-T base pairs), and 
the overstretching force will be where the average free- 
energy change for forming a base-pair is zero. However, 
because of the "average-base" nature of our model, the 
free-energy is roughly linear in the number of base pairs 
in the duplex. When the slope is positive, the fully asso- 
ciated duplex is most stable, whereas when the slope is 
negative the fully unpeeled state is most stable. Thus, by 
estimating the force for which the slope is zero, we obtain 
a value of F OVCI ps 77 pN for our model at 23° C. This value 
is satisfyingly close to the oft-quoted 65 pN for room 
temperature overstretching, and much closer than has 
been previously achieved for any coarse-grained model 
or atomistic simulation. More precisely, at 500 mM salt, 
the concentration at which our coarse-grained model has 
been fitted, Wenner et at. obtain a value of 67pNP and 
and in the more recent study of Zhang et at. the onset 
and mid-point of the transitions were found to occur at 
65 pN and 68 pN, respectively.^ 

The same approach can be used to study the temper- 
ature dependence of F ovor . As the temperature is in- 
creased, the duplex is destabilized relative to ssDNA, so 
i^ovcr is expected to decrease. This behaviour is exactly 
what we see for our model in Fig. [7] Also included in 
this figure is the experimental data of Refs. [7] and [T3] 
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FIG. 6. Free energy profiles as a function of the number 
of base pairs for different forces at (a) room temperature 
(T = 23°) and (b) T = 94°C. The latter is just above the 
zero-force bulk melting temperature and exhibits reentrance. 
Note the relative flatness of the free energy landscapes at this 
temperature. 

Our results show a very similar temperature dependence 
over the temperature range studied in experiments, al- 
beit with our results being displaced to slightly higher 
force by a similar amount to that already noted at room 
temperature. In particular, Zhang et at. find a slope of 
— 0.44pNK _1 at and above room temperature^ which 
is similar to the — 0.48pNK _1 that we obtain for our 
model. Interestingly, Zhang et al. find a change in the 
sign of the slope below 20° C which they interpret as a 
crossover to an S-DNA mechanism. Our results show no 
such change in slope, consistent with the fact that over- 
stretching always occurs by unpeeling in our model. 

One particularly interesting feature at high temper- 
ature and low force is the non-monotonic behaviour of 
F OVCI (T). There is a narrow range of temperature just 
above the bulk (i.e. infinitely long strands) zero-force 
melting temperature where there are two transitions as 
a function of force. This behaviour is illustrated in the 
inset to Fig. [7] and also in Fig. [6jb) where there are two 
changes of sign in the slope of the free-energy profile as a 
function of force for a temperature in this window. At low 
and high force it is favourable for the duplex to unpeel, 
but there is an intermediate force range (~2.5-10 pN at 
94°C) where duplex DNA is most stable. This type of 
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FIG. 7. Dependence of the overstretching force on tem- 
perature. As well as the results for our model, we include 
experimental results from Refs. [7| and [13] at a salt concen- 
tration of 500 mM. Note that the results of Zhang et al. are 
for the onset (not the midpoint) of the transition. The inset 
provides an expansion of the high temperature region where 
non-monotonic behaviour is observed. The horizontal line in 
the inset is the force at which ssDNA and dsDNA have the 
same extension. 



behaviour has been predict ed for a number of theoreti- 
cal modela^ 21 13 1 32 ' 33 1 37 1 38 1 52 1 and occurs because which of 
dsDNA or ssDNA is more extensible depends on the mag- 
nitude of force. At low force dsDNA is more extensible 
because it has a larger persistence length and so has a 
smaller entropy cost for alligning with the force. How- 
ever, when the extension approaches the contour length 
of dsDNA the force required to further extend the duplex 
grows rapidly, and ssDNA becomes more extensible. The 
crossover between these two regimes occurs at approxi- 
mately 6.5 pN (i.e. when the two force-extension curves 
in Fig. [3][a) cross). For forces below 6.5 pN, dsDNA is 
stabilized with respect to ssDNA by increasing force, but 
above this value of the force, the opposite is true. Hence, 
the turning point in F ovm (T) is expected to occur at this 
force. 

One of the features in experiments when overstretch- 
ing occurs by unpccling is that the transition is not fully 
reversible, and that on reducing the force from above 
Fever the system does not trace out exactly the same 
force-extension curve. One can envisage two possible 
sources for this lack of reversibility. Firstly, if the un- 
peeled strand has formed some secondary structure (e.g. 
hairpins), the intramolecular base pairs involved would 
first have to be broken. Secondly, even for an unstruc- 
tured unpeeled strand the reassociation dynamics have 
the potential to be slow, particularly if the strands are 
long, as is typically the case in experiments. 

In Fig. |8ja) we show force-extension curves when the 
force is first increased and then decreased after reach- 
ing different stages in the overstretching process. If the 
force is decreased after overstretching when the non- 
force-bearing strand has dissociated, the resulting force- 
extension curve just follows that for ssDNA (Case A), 
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FIG. 8. (a) Force-extension curves illustrating the irre- 
versibility in our model. In all the simulations the magni- 
tude of the pulling rate is 2.43 x 10~ 6 pN/time step. The 
increasing force run starts from a fully-bonded duplex, and 
provides the starting configurations for the decreasing force 
simulations. In case A, the force is decreased after reaching 
146 pN, whereas in cases B and C the force is decreased after 
reaching 85.3 pN. In case C intramolecular base pairing is not 
allowed, preventing secondary structure formation in the un- 
peeled chain, (b) The final zero-force configuration for case 
B. 



because reassociation of the separated chains does not 
occur. Also depicted in Fig. are two cases when 
the force decrease starts roughly halfway through the 
overstretching transition. In both cases the extension, 
perhaps somewhat surprisingly, initially continues to in- 
crease. As mentioned earlier, the dynamic pulling sim- 
ulations somewhat overshoot the overstretching transi- 
tion. So, even though the force is decreasing, initially 
it is still above the equilibrium overstretching force and 
so the molecule continues to unpeel. However, although 
the extension nearly reaches that for ssDNA, the chains 
do not dissociate, because, as noted earlier, the last few 
base pairs are harder to melt. As the force decreases fur- 
ther, the two cases show quite different behaviour. In 
case B reassociation does occur but only partially and in 
"bursts" . Fig. [^b) depicts the final configuration that 
results, clearly showing that full reassociation has been 
blocked by secondary structure formation in the non- 
force-bearing strand. Similarly, the bursts occur when 
sections involved in secondary structure open up, allow- 
ing further reassociation. 

By contrast, in case C we prevent secondary struc- 
ture in the unpeeled strand by turning off intramolec- 
ular base pairing. Now when reassociation occurs, it 
does so relatively quickly and goes to completion, al- 
beit at a lower force than for overstretching. Thus, for 
our model secondary structure formation is a much more 
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FIG. 9. (a) Force-extension curve for pulling scheme II ob- 
tained from Monte Carlos simulations compared to that for 
scheme I. There are no points above the midpoint of over- 
stretching because the system can no longer bear a force after 
dissociation, (b) Snapshot near to the middle of the transition 
at F= 78 pN. 



substantial contributor to irreversibility than the under- 
lying dynamics of reassociation. This conclusion is an 
agreement with Gross et al. who found that prolonged 
stalls in the reannealing occurred at points at which par- 
ticularly stable secondary structure was possible in the 
unpeeled strand. 27 However, we should note that our sys- 
tem is much smaller than those typically studied in over- 
stretching experiments (for example, A-DNA, which is 
often used in experiments, has 48 502 base pairs). 

Having considered pulling scheme I in considerable de- 
tail, we also examine the other two pulling schemes in 
Fig. [2] Firstly, for scheme II we expect a similar be- 
haviour to scheme I — there are still two free ends that 
allow unpeeling — except beyond the overstretching tran- 
sition, where after dissociation of the two strands the sys- 
tem will no longer be able to bear a force. In Fig. [9] we 
see that the overstretching transition occurs at approxi- 
mately the same force as for scheme I, and again occurs 
by unpeeling. Note that there are no points above the 
midpoint of overstretching due to this instability. 

Scheme III is more interesting, because now over- 
stretching can no longer occur by unpeeling. Two re- 
maining possible mechanisms for overstretching are force- 
induced melting by bubble formation, and an S-DNA-like 
transition to an overstretched but still associated du- 
plex form. The former is expected to require a larger 
force than unpeeling because both strands are still force- 
bearing after bubble formation and so have less entropy 
than an unpeeled strand. Interestingly, in room temper- 
ature experiments at 150 mM NaCl where unpeeling can- 
not occur and the system is not torsionally constrained, 
reversible overstretching still occurs at about 65 pN.^ 
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FIG. 10. (a) Force-extension curves and (b) snapshots of 
the DNA configuration at F = 116 pN as it undergoes over- 
stretching for pulling scheme III obtained from Monte Carlo 
simulations. For comparison, curves for pulling ssDNA and 
dsDNA in pulling scheme I have been added in (a). 



Fig. |10[a) shows the force-extension curve for our 
model for pulling scheme III. The clear difference from 
scheme I and II is that overstretching now occurs at sig- 
nificantly higher force, namely at about 110-115 pN. Al- 
though overstretching can no longer occur by unpeeling, 
Fig. 10]jb) shows that it is still by force-induced melt- 
ing, but in this case by bubble formation. Notably, even 
though the duplex reaches higher forces than in scheme I 
and II, there is still no sign of a transition to an S-DNA- 
like overstretched form. Also, that the force-extension 
curve above overstretching does not follow that for ss- 
DNA is simply because each individual strand only bears 
a force F/2 (Fig. [2jc)) , and therefore the curve can be 
made to overlap with that for ssDNA by shifting it down 
by a factor of 2. 



IV. CONCLUSIONS 

In this paper we study the overstretching transition 
in a recently developed coarse-grained model of DNA, 
which has previously been shown to give an excellent de- 
scription of the thermo dynam ic, structural and mechan- 
ical properties of DNA.EM1 At room temperature and 
500 mM salt our model undergoes overstretching by un- 
peeling at 77 pN, just somewhat higher than the experi- 
mentally observed overstretching force, and much closer 
than the predictions for any other coarse-grained poten- 
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tial modeP^Hlor any atomistic simulations.^HSU We are 
also able to reproduce the temperature dependence of 
the overstretching force well. This ability to accurately 
reproduce phenomena to which it was not directly fitted 
provides further validation of our model. 

It is also important that we understand why our coarse- 
grained model does so well in describing the overstretch- 
ing transition. As emphasised in Section |II C| to predict 
force-induced melting correctly, one just needs to know 
the zero-force thermodynamics of melting and the force- 
extension curves of dsDNA and ssDNA. As our model 
describes all of the above reasonably well, its success in 
describing overstretching by force- induced melting is per- 
haps not surprising. 

The 9-10 pN overestimation of the room temperature 
overstretching force in our model is probably due to a 
number of factors. Firstly, our description of the ss- 
DNA force-extension curve deviates slightly from fits to 
experiment,^ in particular, underestimating the exten- 
sion at larger forces. The net effect is to somewhat un- 
derestimate the stabilization of ssDNA by force near to 
overstretching. Secondly, the zero-force thermodynamics 
in our model has been designed to fit the thermodynamic 
predictions of the SantaLucia nearest-neighbour modelj^ 
and so any errors in this model will be carried over to 
our results. For example, Huguet et al. have shown 
that using the SantaLucia model and fits to the force- 
extension curves of ssDNA and dsDNA gives an approx- 
imately 1 pN overestimation of DNA unzipping at room 
temperature.^ As the difference in extension between 
the duplex and ssDNA is less for the overstretching than 
for the unzipping geometry, one would expect the error 
due to the use of the SantaLucia model to be larger for 
overstretching, with the above result implying an over- 
estimation of the order of 4pN. However, when Gross et 
al. applied a similar approach to model their experimen- 
tal overstretching force-extension curves, the model actu- 
ally underestimated the experimental results.^ Finally, 
in order to make the computation of the overstretch- 
ing force tractable, we turned off intramolecular base- 
pairing. This base-pairing stabilizes the ssDNA state and 
so any effect of this approximation will be to contribute 
to an overestimation. 

In all simulations where free ends are present, we see a 
force-induced melting mechanism for overstretching oc- 
curring by unpeeling, one of the confirmed experimen- 
tal mechanisms.^ However, no sign is ever seen of an 
S-DNA-like transition to an overstretched duplex form. 
This shortcoming leads to an overstretching force when 
there are no free ends (i.e. scheme III) that is much higher 
than experiment, because the mechanism is force-induced 
melting by bubble formation. The reason for our model's 
inability to exhibit an S-DNA-like state is harder to iden- 
tify, as not much is known about the detailed structure 
of S-DNA. One suggestion is that S-DNA might adopt an 
unstacked but base-paired configuration!^^ However, if 
this were the case one might expect our model to be able 
to reproduce this, as it describes the thermodynamics of 



stacking and base-pairing well. Furthermore, the reason 
why such a form is not stable in our model is simply that 
stacking provides a larger free energetic contribution to 
the stability of B-DNA than base pairing and that the 
overstretched unpeeled state retains a significant frac- 
tion of its stacking (Fig. [5| . This stacking occurs both 
in the non-force-bearing unpeeled strand that is totally 
free to stack and in the force-bearing strand where short 
runs of stacked bases can align their backbones along the 
force axis. Even when there are no free ends, for which 
force-induced melting in our model occurs at higher force 
by bubble formation with less stacking retained because 
both strands are force-bearing, an unstacked but base- 
paired duplex is still not observed. Our results strongly 
suggest that S-DNA does not correspond to such a struc- 
ture, but instead probably to a more exotic structure. A 
plausible candidate structure for S-DNA might be "zip- 
DNA" j 43 l 49 l where the chains are no longer base-paired 
but interdigitate allowing stacking between bases on dif- 
ferent strands. Our model does not include the possibility 
of this more unusual stacking interaction, and so of course 
would not be able to reproduce any such behaviour. 
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